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Abstract 

We consider a general stochastic branching process, which is relevant to earthquakes, and study the 
distributions of global lifetimes of the branching processes. In the earthquake context, this amounts to the 
distribution of the total durations of aftershock sequences including aftershocks of arbitrary generation 
numbers. Our results extend previous results on the distribution of the total number of offsprings (direct 
and indirect aftershocks in seismicity) and of the total number of generations before extinction. We 
consider a branching model of triggered seismicity, the ETAS (epidemic-type aftershock sequence) model 
which assumes that each earthquake can trigger other earthquakes ( "aftershocks" ) . An aftershock sequence 
results in this model from the cascade of aftershocks of each past earthquake. Due to the large fluctuations 
of the number of aftershocks triggered directly by any earthquake ("productivity" or "fertility"), there 
is a large variability of the total number of aftershocks from one sequence to another, for the same 
mainshock magnitude. We study the regime where the distribution of fertilities [i is characterized by a 
power law ~ 1 //x 1+7 and the bare Omori law for the memory of previous triggering mothers decays slowly 
as ~ l/t 1+s , with < 6 < 1 relevant for earthquakes. Using the tool of generating probability functions 
and a quasistatic approximation which is shown to be exact asymptotically for large durations, we show 
that the density distribution of total aftershock lifetimes scales as ~ l/t 1+9 ^ when the average branching 
ratio is critical (n = 1). The coefficient 1 < 7 = b/a < 2 quantifies the interplay between the exponent 
b ~ 1 of the Gutenberg-Richter magnitude distribution ~ 10~ bm and the increase ~ \0 arn of the number 
of aftershocks with the mainshock magnitude m (productivity) with a ~ 0.8. The renormalization of the 
bare Omori decay law ~ l/t 1+e into ~ l/t l+e ^ stems from the nonlinear amplification due to the heavy- 
tailed distribution of fertilities and the critical nature of the branching cascade process. In the subcritical 
case n < 1, the cross-over for ~ l/t l+e ^ at early times to ~ l/t l+e at longer times is described. More 
generally, our results apply to any stochastic branching process with a power-law distribution of offsprings 
per mother and a long memory. 



* Electronic address: sornette@moho.ess.ucla.edu 
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I. INTRODUCTION 



We study the distribution of the total duration of an aftershock sequence, for a class of branch- 
ing processes appropriate in particular for modeling earthquake aftershock sequences. The 
noteworthy particularity and challenging property of this class of branching processes is that the 
variance of the number of progenies in direct lineage from the mother is mathematically infi- 
nite. In addition, a long-time (power law) memory of the impact of a mother on triggering her 
first-generation daughters gives rise to subdiffusion jjj |^ and non-mean field behavior in the dis- 
tributions of the total number of aftershocks per mainshock and of the total number of generations 

n 

before extinctions |5|]. Here, we add on these previous works but showing that the distribution of 
the total duration of an aftershock sequence is extremely long-tailed: the very heavy-tailed nature 
of the distribution of the durations of aftershock sequences predicted by this simple model may 
explain the large variability of the lifetimes of observed aftershock sequences and is compatible 
with the observation that felt aftershocks of the great Mino-Owari (1891) Japanese earthquake, 
that inspired Omori's statistical rate model, have persisted at a rate consistent with the Omori law 
for 100 years [6|. 

Our results may also be of interest to other systems which are characterized by branching pro- 
cesses with a broad power-law distribution of fertilities, such as epidemic transmission of diseases, 
and more generally transmission processes involving avalanches spreading on networks such as the 
World Wide Web, cellular metabolic network, ecological food webs, social networks and so on, 
as a consequence of the well-documented power law distribution of connectivities among nodes. 
Our results are thus relevant to systems in which the number of offsprings may be large due to 
long-range interactions, long-memory effects or large deviation processes. 



II. THE EPIDEMIC-TYPE AFTERSHOCK SEQUENCE (ETAS) BRANCHING 
MODEL OF EARTHQUAKES WITH LONG MEMORY 

We consider a general branching process in which each progenitor or mother (mainshock) is 
characterized by its conditional average number 

N m = «//(m) (1) 
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of children (triggered events or aftershocks of first generation), where 

/i(m) = I0 a ( m ~ m o) ; (2) 

is a mark associated with an earthquake of magnitude m > mo (in the language of "marked point 
processes"), k is a constant factor and mo is the minimum magnitude of earthquakes capable 
of triggering other earthquakes. The meaning of the term "conditional average" for N m is the 
following: for a given earthquake of magnitude m and therefore of mark p(m), the number r of its 
daughters of first generation are drawn at random according to the Poissonian statistics 

^(r) = ^e-^ = ^e--. (3) 

Thus, N m is the expectation of the number of daughters of first generation, conditioned on a 
fixed magnitude m and mark p(m). The expression (j2J) for p,(m) is chosen in such a way that it 
reproduces the empirical dependence of the average number of aftershocks triggered directly by 
an earthquake of magnitude m (see []| and references therein). Expression (JTJ) with (J2J) gives the 
so-called productivity law of a given mother as a function of its magnitude. 

In addition, we use the well-known Gutenberg-Richter (GR) density distribution of earthquake 
magnitudes 

p{m) = b ln(10) l()- b ( m - m °) , m > m , (4) 

such that p(x)dx gives the probability that an earthquake has a magnitude equal to or larger 
than m. This magnitude distribution p(m) is assumed to be independent on the magnitude of the 
triggering earthquake, i.e., a large earthquake can be triggered by a smaller one [?], 

Combining (jlj) and (J2J), we see that the earthquake marks p and therefore the conditional average 
number N m of daughters of first generation are distributed according to a power law 

PM = ^x&> 1<//<+oo, 7 = 6/a. (5) 

Note that p^p) is normalized: J^°° dp, p^p) = 1. For earthquakes, 6 ~ 1 almost universally and 
a ~ 0.8 |2(, giving 7 w 1.25. The fact that 1 < 7 < 2 implies that the mathematical expectation of 
p and therefore of N m (performed over all possible magnitudes) is finite but its variance is infinite. 

For a fixed 7, the coefficient k then controls the value of the average number n of children of 
first generation per mother: 

n = (N m ) = k(h) = k— — - , (6) 

7-1 
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where the average (N m ) is taken over all mothers' magnitudes drawn from the GR law. In the 
terminology of branching processes, n is called the branching ratio. For n < 1, there are on average 
less than one child per mother: this corresponds to transient (sub-critical) branching processes with 
finite lifetimes with probability one. For n > 1, there are more than one child per mother: this 
corresponds to explosive (super-critical) branching processes with a number of events growing 
exponentially with time. The value n — 1 of exactly one child per mother on average is the critical 
point separating the two regimes. 

Finally, we assume that a given event (the "mother" ) of magnitude m > m occurring at time 
ti gives birth to other events ("daughters") of first generation in the time interval between t and 
t + dt at the rate 

Mt) = N m Ht - U) = N m Hit) (7) 

where < 9 < 1, H(t) is the Heaviside function, c is a regularizing time scale that ensures that 
the seismicity rate remains finite close to the mainshock and N m is given by The time decay 
rate is called the "direct Omori law" [ill Q • -Due to the process of cascades of triggering by 
which a mother triggers daughters which then trigger their own daughters and so on, the direct 
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which is the 



Omori law (J7j) is renormalized into a "dressed" or "renormalized" Omori law 
one observed empirically. 

Expressions (I1I2I4I7J) define the Epidemic-Type Aftershock Sequence model of triggered seismic- 
ity introduced by Ogata in the present form |9| and by Kagan and Knopoff in a slightly different 
form 11 01. 



III. GENERAL FORMALISM IN TERM OF GENERATING FUNCTIONS 

Since we are interested in characterizing the distribution of the random times at which an 
aftershock sequence triggered by a given mainshock terminates, we take the time of the mainshock 
of magnitude m at the origin t = and we do not consider the effect of earlier earthquakes. 
This is warranted by the fact that sequences of earthquakes generated by different mainshocks are 
independent in the ETAS branching model. 
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A. First generation aftershocks 

Let us first discuss some more detailed statistical description of first generation aftershocks. 
Each aftershock arising independently from another preceding aftershock itself born at the random 
time ti has its birth time possessing the probability density function (PDF) $(t — £j) defined in 
(J7J) and cumulative distribution function (CDF) b{t) = J^§{t')dt'. Here and everywhere below, 
the dimensionless time t/c is used and we replace t by t/c, with the understanding that t or r 
means t/c when needed. It is convenient to introduce the complementary CDF of first generation 
aftershocks 

a{t) = 1 - b(t) = . (8) 

Let us consider a mainshock with mark /x that triggers exactly r aftershocks of first-generation 
arising at the moments (ti,t 2 , . . . ,t r ). Then, the CDF of the time T(fi\r) of the last arising 
aftershock is equal to 

P^(t\r) = Pr (T(/x|r) = max^, t 2 , . . . , t r } < t) = [b(t)] r . (9) 

Averaging this CDF over the random first-generation aftershock numbers r at fixed /i weighted 
by their probability p^(r) given by © yields the CDF P^(t) for the total duration T(fi) of the 
first-generation aftershocks 

P„(t) = Pr (T(ai) < t) = G${t)] . (10) 

Here, G^{z) = J2'^LoP^( r ) zr * s the generating probability function (GPF) of the number of first- 
generation aftershocks. For the Poissonian statistics (J3j), it is equal to 

G,(z) = e K ^ z - 1 K (11) 

This leads to the well-known relation 

P M (t) = <,-«/"»(*) . (12) 

In the ETAS model, the Gutenberg-Richter distribution (@J) of magnitudes together with the 
productivity law (j2J) implies the power law (jHJ) for the marks fi. Averaging over all possible main- 
shock magnitude thus amounts to averaging (fTUj) over all possible /i's. The CDF of durations T of 
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first-generation aftershocks generated by some mother of arbitrary magnitude arising at time t — 
is equal to 

P (t) = G[b(t)} , (13) 

where G(z) = (G^z)) is the average of G^[b(t)} over the random magnitudes m (or equivalently 
random marks /i) In the relevant case of the Poissonian GPF (jllj) and using (j3J), we obtain 

G(z)= 1 ^(l-zrT(- 7} K(l-z)) } (14) 

where T(x,y) is the incomplete Gamma function and 7 = b/a. For real aftershocks, 1 < 7 < 2 
and a typical value is 7 ~ 1.25. Then, it is easy to show that the first terms of G(z) in a power 
expansion with respect to 1 — z are 

G(z) ~ 1 - n(l - z) + (3(1 - zf , 1<7<2, (15) 

with n given by © and 

p = li ,(2zlVlQz±. (16) 

V 7 / 7-1 

B. All generation aftershocks 

In the ETAS model, any event (the initial mother or any aftershock, whatever its generation 
number) triggers its aftershocks of first-generation in a statistically independent and equivalent 
manner, according to the laws given in section [H] This gives the possibility of obtaining closed 
equations for the CDF of the total duration of aftershocks triggering processes. 

Let T be the random waiting time between a mainshock and one of his first-generation af- 
tershocks, chosen arbitrarily. The PDF of T is nothing but $(t) defined in ((7|). Let T be the 
random duration of the aftershock branching process triggered by this first-generation aftershock. 
The CDF of T is denoted P(t). Then, the total duration, measured since the mainshock, of the 
sequence of aftershocks generated by this pointed out first-generation aftershock is T + T. The 
CDF ¥(t) of this sum is therefore the convolution 

F(t) = ® P(t) . (17) 

Replacing in (fTUj) b(t) by F(i) and taking into account the equality (fTTj) . we obtain the CDF of 
the total duration T(/z) of a sequence of aftershocks over all generations of a given event of mark 



\i that occurred at t — 0: 

P M (t) = Pr (T(/i) < t) = e -^ MW , (18) 

where 

R(t) = 1 - F(t) (19) 

is the complementary to the F(i) CDF defined in ()17|) . Correspondingly, replacing in (|13j) P (t) by 
P(t) and &(£) by F(£), we obtain the self-consistent equation for the CDF ¥{t) 

P(t) = G[F(t)] = G [$(t) ® P(t)] . (20) 

It is convenient to rewrite (J2(Jj) as 

R(t) - Q(t) = n [R(t)} , (21) 

where Q(t) = 1 - P(i) and 

= G(l-z) +2- 1. (22) 

For our subsequent analysis, expression (|2T|) is more convenient than equation (j271|l for the following 
reasons. First of all, instead of the CDF's P(i) and F(i) entering in (|2())1. equation (|2*T|) is expressed 
in terms of the complementary CDF's Q(t) and R(t) which both tend to zero for t — > oo. In 
addition, the function Q(z) also tends to zero for z — > 0. This gives the possibility of extracting 
the influence of the nonlinear terms of the GPF G(z) on the asymptotic behavior of the solution 
for t — > oo. Indeed, at least for 7 < 1.5, the GPF G(z) is very precisely described by the truncated 
series (115)1 . The corresponding series for Q(z) is 

n(z) ~ (1 -n)z + /3z 7 , (23) 

which reduces to a pure power law in the critical case n = 1: 

Q{z) ~ /52 7 . (24) 

Correspondingly, in the critical case n = 1 and most important for earthquake applications for 
which 1 < 7 < 2 holds, equation (|2*T|) has the form 

R(t) - Q(t) = (3R y (t) . (25) 

The exact auxiliary function Q(z) defined by ()2*2*j) for n = 1 and its power approximation for 
7 = 1.25 are shown in Fig. 1. 

Our goal is now to solve (|2*T|) and in particular (|25j) to explore in details the statistical properties 
of the durations of aftershocks sequences, resulting from cascades of triggered events. 
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IV. FRACTIONAL ORDER DIFFERENTIAL EQUATION FOR THE COMPLEMEN- 
TARY CDF R(t) 



In order to exploit equation (J21|) . we first need to express Q(t) as a function of For this, 

we note that expression (|T7|) is equivalent to 

R(t) = a(t) + ® Q(t) , (26) 

as can be seen from direct substitutions using (|SJ), (|T9"j) and Q(i) = 1 — P(£). Applying the Laplace 
transform to both sides of this equality, one gets 

a, , R(s) l-$(s) 



$(s) s$(s) 

where 



(27) 



/•OO 

$( s )=/ $(t)e- s *rft = ^(cs) 9 e cs r(-^,s), (28) 
Jo 

where we have made the correspondence t — > t/c explicit (where c is defined in (JJJ)). We shall 
be interested in the probability distribution of the durations of total sequences of aftershocks for 
durations much larger than c. In this case, one can replace $(s) by its asymptotics for small s 

^s)^l-S(cs) 9 ^ 1 + ^ cs)e , ca<l, (29) 

where 5 = T(l — 6). Substituting it into (J27|) leads to 

Q(s) = [1 + 5(c s) 9 ] R(s) - 6 cV" 1 , (30) 

which is equivalent, under the inverse Laplace transform, to the fractional order differential equation 

, m+s ^-n\ ( 31 ) 



dt e Vt, 

Equation (|21j) thus yields the following fractional order differential equation for K.(i) (going back 
to the reduced time variable r = t/c) 

5^ + fl(R) = r' e . (32) 
In particular in the critical case n — 1, using the power approximation ([24)1 . we obtain 

50 + / ^ = r~ e . (33) 
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Note that the nonlinear fractional order differential equation (|32|) is exact for $(t) given by 

M) = -L^e (^js) , (34) 



where &g(t) is the fractional exponential distribution possessing the following Laplace transform 

= rib*' (35) 

which has the integral representation 

POO -t , . 

$e(r) = / -exp( — ) £ e (x)cte, (36) 
Jq x \ xJ 

where 

* / s 1 sin(7r6>) .„„. 

TlX X" + X + 2 COS{ir9) 

One can interpret (|3fi|) as the decomposition of the fractional exponential law into regular expo- 
nential distributions, and £e(x) given by (|37j) as the "spectrum" of their mean characteristic decay 
time x. For 9 — > 1, the spectrum (|3*7|) weakly converges to the delta-function 5(x — 1) and the 
fractional exponential law transforms into the regular exponential distribution $i(r) = e~ T . For 
9 = 1/2, there is an explicit expression for the fractional exponential distribution 



$1/2 (t) = \ — - e T erfc (y/r) . 

V TIT 



(38) 



It is easy to show that the asymptotics of the fractional exponential distribution are 



r 9 - 1 9T- - 1 



W-fTm $e(r) ^ r(l-g) (T>>1) - (39) 

Fig. 2 shows a log- log plot of the Omori law $(£) defined in (jZJ and of the corresponding 
fractional exponential distribution (|3*4*j) as a function of the reduced time r = t/c and for # = 1/2, 
demonstrating the closeness of these two distributions. 



V. EXACTLY SOLUBLE CASE: PURE EXPONENTIAL OMORI LAW 

Before addressing the case of interest of earthquakes where the direct Omori law $(t) is a power 
law with exponent < 9 < 1, it is instructive to present the solution for the case where $(£) is an 
exponential. In this case, an exact solution can be obtained in close form. This exact solution will 
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be useful to check the quasistatic and dynamical linearization approximations developed below to 
solve the difficult case where $(£) is a power law with exponent < 6 < 1. 
We write the exponential direct Omori law in non-reduced time as 

*(t) = ±e*p(--) =► *(*) = —!-, (40) 
c \ c J I + cs 

so that equation (|2*7)l transforms into 

Q(s) = (l + cs)E(s) -c. (41) 
After inverse Laplace transform, we get 

Q(t) = R(t) + - cS(t) , (42) 



and equation (|2*T|) takes the form 

dR(t) 



n [M(t)] = , (43) 



dt 

or, in the more traditional form of a Cauchy problem 

^ + O [E] = , E(r = 0) = 1 . (44) 

GST 

The numerical solution of (jH| is easy to obtain. In addition, using for Q(z) the series approximation 
(|23p. one obtains the analytical solution of the Cauchy problem f!44|) under the form 



1 + exp (1 - n)- 



1 — n J \ g J 1 — n 

where g — 1/(7— 1). In particular, in the critical case n — 1, this leads to 



(45) 



1 + ^r^l \ (46) 





Fig. 3 shows the numerical solution of equation (|44j) together with its analytical solution (|45|) 
obtained using the polynomial approximation (J23)) of the function Q(z) defined in (|22j) . for 7 = 1.25 
and n = 0.99. It is seen that curves are very close each other. 

Note that, in the sub critical case n < 1, there is a crossover from the power law ()46|) at early 
times which is characteristic of the critical regime n — 1, to an exponential decay at long times of 
the complementary CDF R. 
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VI. DYNAMICAL LINEARIZATION AND QUASISTATIC APPROXIMATIONS TO 
OBTAIN THE ASYMPTOTIC TAIL OF THE DISTRIBUTION OF TOTAL AFTER- 
SHOCK DURATIONS 



A. Linear approximation 

To obtain some rough estimate of the complementary CDF R(i), let us consider the linearized 
version of the fractional order differential equation (|32|) 

5^ + v R = t-\ (47) 

where the following linearization has been used 

Q[R] ~r]R, r] = £1(1) = G(0) . (48) 

The Laplace transform of the solution of the linearized equation (J47j) has the form 

Ms) = -^TT ■ (49) 
r] + os 

The corresponding complementary CDF is equal to 

R = E (-^r ), 6 = T(l-e), (50) 



s 

where Eg(z) is the Mittag-Leffler function. Its integral representation is 

f-oo y e - 1 e y dv 

E e (-x) = -smne — \ U a (x > 0) . (51) 

In particular for 9 = 1/2, it is equal to 

E 1/2 {-x) = e^erfc (x) . (52) 

Its asymptotics reads 

EJ-T.) ~ 

xS 



E e (—x) ~ — (a; — > oo) , (53) 



which is already very precise for x > 2. 

The suggested dynamical linearization approach consists in replacing the factor rj in by 

»« = ^. (54) 
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to correct for the nonlinear decay of the relaxation of the complementary CDF R as a function of 
time. It is interesting to check the validity of this dynamical linearization procedure for the exactly 
solvable exponential Omori law (J4(Jj) . In this case, the solution of the linearized equation (J44j) is 



-r/r 



(55) 



Substituting here (|54|) for 77, we obtain in the critical case the following transcendent equation 



exp(— t (3 ' 



(56) 



Its solution is equal to 



where 



Y(x) 



x 



1 



9 



x 



rl3 



(57) 



(5f 



7-1 9 

and Y(x) is the solution of the transcendent equation Y e Y = x. For x > 2, there is very precise 
approximate solution of this equation: 



Y(x) ~ lnx 



l + (l + lnx) 1-W1 + 



2 ln(lnx) 
(1 + lnx) 2 



lnx . 



(59) 



Thus, for large x, the main asymptotics of the dynamical linearization approximation (|57|) of the 
Cauchy problem (|4~4"j) differs from the main asymptotics R ~ x~ 9 of the exact solution (j4T)j) only 
by logarithmic correction In 9 x. 



B. Quasistatic approximation 

Close inspection of the complementary CDF (j5Uj) and its asymptotics 

1 /M\ 1/61 
^ — — j ^>r*, r*= - (60) 

derived from relation (|53|) gives us a hint of how to approach the solution of the nonlinear fractional 
order differential equation (|32|) and (JHHj) by using a quasistatic approximation. Indeed, notice that 
the asymptotics (joT7|) is solution of the truncated equation (j4"7j) 

t)R = t~ 9 , (61) 
13 



where we omitted the fractional order derivative term. 

Applying this same quasistatic approximation to the nonlinear fractional order differential equa- 
tion gives the approximate equality 

~ T- e . (62) 
In particular, in the critical case n — 1 for which Q(z) ~ /3z 7 , we have /3R 7 ~ r~ e , or equivalently 

Expression (|63|) will lead to your main result (|68|) below. 

The validity of this quasistatic approximation is checked by calculating the derivation of frac- 
tional order 9 of the approximate solution (JBlIj) . Using the standard tabulated formula of fractional 
order analysis 

gg = rq+p) ^ 

we obtain 

^-r^BU-^^, (65) 



(ir 61 7 + 1 \ ' 7, 

where B(x, y) is the Beta function. For any fixed 1 < 7 < 2 and < 9 < 1, there is a r*(7, 0) < 00 

such that 

5^ </?M 7 ~r- e if r>r*( 7 ,^) (66) 

so that the quasistatic approximation becomes applicable. The physical background of the power 
asymptotics (fHUj) of the solution of the linear equation (|37j) and of the quasistatic approximation (jfilfy 
of the nonlinear equation is obvious: the asymptotics R ~ r -9 given by ()6())1 is a consequence of 
the power tail $(£) ~ t -6 *" 1 of the bare Omori law, while the more slowly decaying R ~ r^ 9 ^ given 
by ()63|) is the result of an interplay between the long memory property of the bare Omori law and the 
amplification by the power law Q(z) ~ z 1 , a signature of the broad distribution of productivities of 
daughter aftershocks from mother earthquakes. This gives rise to a renormalization of the exponent 
9 into a smaller exponent 9/j (for 1 < 7 < 2). 



C. PDF of the total duration of aftershocks branching processes 



The previous sections have discussed in details how to obtain the complementary CDF R of 
the total duration of aftershock branching processes, corresponding to some first generation after- 
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shock, triggered by a main earthquake. The CDF P M of the total duration of aftershock triggering 
processes, taking into account all aftershocks triggered by a main earthquake of fixed magnitude, is 
described by relation (|18|) . The corresponding PDF of the total duration of an aftershock sequence 
is thus equal to 

w M(r) = _^ e -^w^). (67) 

If /i/t ^> 1 (as is the case for a large earthquake which has a large average productivity), then, due 
to the exponential factor in (JBTj) . this PDF differs significantly from zero only if R is very small. 
Then using the expression for small values of K described by the quasistatic approximation (}63|) . 
we obtain 



Expression (joH|) is our main result. Fig. 4 shows a log- log plot of the PDF (j68|) for different values 
of the mainshock size fiK for 7 = 1.25 and 9 = 0.2 (Recall that /3 is given by ()16j) and we put it 
equal to 1 to draw Fig. 4). 

Expression (}68|) shows that the power law tail holds for durations t/c > t^/c oc (fin) 1 ^ 9 ~ 
X0( Q7 / e ) m for which the exponential factor goes to 1. Thus, for 9 small (~ 0.1 — 0.3 as seem to 
be relevant for earthquakes), expression (|F)K|) exhibits a very strong dependence on the mainshock 
magnitude through its impact (0) on the mark /i. Therefore, the most relevant part of the distribu- 
tion of the durations for small mainshocks is controlled by the power law tail r~ l ~ e ^ . In contrast, 
the observable part of the distribution of durations for very large mainshocks is controlled by the 
exponential term which, together with the power law prefactor, leads to a maximum: for very large 
fi, W At (r) starts from zero for r = and then increase up to a maximum before crossing over slowly 
to the power law tail r- 1 -^, as illustrated in Fie;. 4. 



D. Crossover from critical to subcritical regime 

The asymptotics of the complementary CDF M satisfies the equation 1)62)1 in the quasistatic 
approximation. In the subcritical regime, using the polynomial approximation ()23)1. one can rewrite 
equation ()6"2"|) in the form 

(1 - n)R + f3W = T- e . (69) 
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It is seen from this equality that if R > M c , where 




1 — n 



) 



(70) 



then one can neglect the linear term in the left-hand-side of equality (joT?|) and obtain the power 
law ()63|). typical of the critical regime n — 1. In contrast, if M. < M. c , then the sub critical scenario 
of the complementary CDF R dominates and equality (J69)) gives the subcritical power law 

K~ . (71) 

1 — n 

It follows from (jHEJ) and ()7())1 that the time of the crossover from the critical to the subcritical 
regime is equal to 
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Fig. 1: Plots of exact Q(z) defined by (|22[) (lower curve) and its pure power approximation 
(1241) (upper curve) for 7 = 1.25 and n = 1. 
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Fig. 2: Loglog plots of the direct Omori law <E>(t) defined in Q (lower curve) and of the 
fractional exponential distribution (upper curve) for 9 = 0.5 and c = 1. 




Fig. 3: Plot of the numerical solution of equation (|44|) for the complementary CDF R of 
the total duration of an aftershock sequence and the corresponding analytical approximate 
expression (|45|) for R for the parameters 7 = 1.25 and n = 0.99. 
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Fig. 4: Log-log plots of the PDF (|68|) of the total aftershock sequence durations for a 
mainshock of mark //, with /ik = 2, 5, 10, 15, for the parameter values 7 = 1.25, 9 = 0.2 and 
n = 1. 
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